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C/3 ■ Abstract 



In this paper we outline an algorithmic approach to compute Puiseux series expansions 
for algebraic surfaces. The series expansions originate at the intersection of the surface with 
as many coordinate planes as the dimension of the surface. Our approach starts with a 
polyhedral method to compute cones of normal vectors to the Newton polytopes of the given 
polynomial system that defines the surface. If as many vectors in the cone as the dimension 
of the surface define an initial form system that has isolated solutions, then those vectors 
are potential tropisms for the initial term of the Puiseux series expansion. Our preliminary 
methods produce exact representations for solution sets of the cyclic n-roots problem, for 
n = to 2 , corresponding to a result of Backclin. 
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■ 1 Introduction 

We presented polyhedral algorithms to develop Puiseux expansions, for plane curves in [2] and 
for space curves in [TJ, based on ideas described in [32]. In this paper we explain a polyhedral 
approach to compute series developments for algebraic sets. Although we use the numerical solver 
of PHCpack [31], one may use any solver for the leading coefficients of the series and obtain a 
purely symbolic method. We implemented our methods using Sage [28] . 

We could reduce the treatment of algebraic sets to the curve case by adding sufficiently many 
hyperplanes in general position to cut out a curve on the set. This approach [26J is not flexible 
enough to exploit permutation symmetry as the added general hyperplanes ignore the symmetric 
structure of the polynomial system. 



*This material is based upon work supported by the National Science Foundation under Grant No. 1115777. 
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Although presently we do not have a fully automatic implementation suitable for benchmark- 
ing on a large class of polynomial systems, we have obtained promising results on the cyclic n-root 
systems: 

X +Xl -\ h X n _l = 

X X 1 + XlX 2 H h X n -2X n -l + X n -lX = 

i = 3, 4, . . . , n - 1 : ]J x fc mod n = 

3=0 k=j 
XqX\X2 ■ ■ ■ X n -l -1 = 0. 

The cyclic n-roots system is a standard benchmark problem in computer algebra, relevant to 
operator algebras. We refer to [30J for recent advances in the classification of complex Hadamard 
matrices. In |12| . the close relationship of ([I]) with some systems occurring in optimal design of 
filter banks is stressed. The numerical factorization of the two dimensional surface of cyclic 9- 
roots into 6 irreducible cubics was reported in [25J. Recent results for the cyclic 12-roots problem 
can be found in [23J. 

Surprisingly, while looking to develop Puiseux series for algebraic sets, for cyclic 9-roots we 
found exact results: the first term of the series satisfies the entire polynomial system. These exact 
result correspond to known (see e.g. [4] or [12]) configurations of cyclic n-roots. 

The type of polynomial systems targeted by the polyhedral approach are sparse polynomial 
systems. We introduce our approach in the next section with a very particular sparse class of 
systems. We use unimodular transformations to work with points at infinity. The second section 
ends with a general approach to solve a binomial system. 

To find the initial coefficients in the Puiseux series we look for initial form systems, systems 
that have fewer monomials than the original systems and that are supported on faces of the 
Newton polytopes. Faces of the Newton polytopes that define the initial forms are determined 
by their inner normals. Those inner normals that define the initial form systems are the leading 
powers (called tropisms) of generalized Puiseux series. The leading coefficients of the series vanish 
at the initial form systems. 

In the third section we define initial form systems, give an illustrative example, and describe 
the degeneration of a (i-dimensional algebraic set along a path towards the intersection with 
the first d coordinate planes. Polyhedral methods give us cones of pretropisms and initial form 
systems that may lead to initial coefficients of Puiseux series. We end this paper giving an exact 
description of positive dimensional sets of cyclic n-roots. 

Related work. A geometric resolution of a polynomial system uses a parameterization of the 
coordinates |14| for global version of Newton's iterator [9]. Our algorithms arose from an un- 
derstanding of [6l Theorem B] and are inspired by tropical methods [7] and in particular by the 
constructive proof of the fundamental theorem of tropical algebraic geometry |20| . Puiseux series 
occur perhaps most often in the resolution of singularities, [3] describes an extension of Newton's 
method using the notion of tropical variety. Software related to [20] is Gfan [18] and the Singular 
library tropical. lib [19] . 

Connections with Grobner bases are described in [29J. Polyhedral and tropical methods 
for finiteness proofs in celestial mechanics are explained in |16] and |17j . Truncations of two 
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dimensional varieties are studied in [21]. The unimodular coordinate transformations are related 
to power transformations in [8]. A Newton-Puiseux algorithm for polynomials in several variables 
is described in [5]. In |22| . fractional power series solutions are developed for generic systems. 
Acknowledgements. We thank Marc Culler for mentioning the Smith normal form for uni- 
modular transformations. We appreciate the comments of the reviewers. 



2 Binomial Systems 

We aim to solve sparse polynomial systems, systems of polynomials with relatively few monomials 
appearing with nonzero coefficient. The sparsest polynomial systems which admit solutions with 
nonzero values for all coordinates consist of exactly two monomials in every equation and we call 
such systems binomial systems. See e.g.: [10] and for more on binomial ideals. 

To represent a d-dimensional solution set S intersecting the first d coordinate planes in as 
many regular isolated points as the degree of S, the first d variables can serve as independent 
parameters. The parameterizations that are of interest to us start with the generators of cones 
of normal vectors defining initial forms of polynomial systems. 



2.1 An Example 

Consider for example 



4^3 



X qX ~^X <^X g 



l 



xqxix 2 xz - 1 

We write the exponent vectors in the matrix 




0. 



A 



2 14 3 
1111 



(2) 



(3) 



and we look for a basis of the null space of A. Two linearly independent vectors that satisfy Ax = 
are for example u = (—3,2,1,0) and v = (—2,1,0,1). Placing u and v in the columns of a 
matrix M leads to a coordinate transformation: 



M 



-3 


-2 


1 





2 


1 





1 


1 














1 









X2 
X 3 



2/0 

yi- 



(4) 



The coordinate transformation x = y eliminates yo and y± — because u and v are in the null 
space of A — as substituting the coordinates corresponds to computing ylu and Av, reducing 
the given system to 

1 = 
1 = 0. 



2/22/3 



(5) 



Solving the reduced system in ([5]) gives values for 2/2 and yz which after substitution in the 
coordinate transformation in (j4|) yields an explicit solution for the original system in ([2]) with yo 
and y± as parameters. 
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2.2 Unimodular Transformations 



In the previous section we constructed in (|4|) a unimodular coordinate transformation x = y M , 
where det(M) = ±1. In the new y coordinates all points that make the same inner product of 
the ith. row of the given exponent matrix A will have the same value for y^. 

The null space of the matrix A is stored in the rows of the matrix B: AB T = 0. The 
Smith normal form of B consists of the triplet (U, S, V), where U and V are unimodular matrices 
(det(U) = ±1 and det(V) = ±1), and the only nonzero elements of S are on the diagonal: 
UBV = S. 

If U equals the identity matrix, then UBV = S implies B = SV~ l . This means that for 
any x, the outcome of Bx. is the same as SV^x. If moreover S contains the identity matrix, 
then V -1 defines the unimodular transformation M. The next examples illustrates the case of 
general U but where S contains the identity matrix. 

For the matrix A in (J3J), the matrix B has in its two rows the vectors u and v so that AB T = 0: 



B 



-3210 
-2101 



(6) 



The computation of the Smith normal form of B with GAP [15J (from the console in Sage [28j 



gives 



and 



U 



V 



s 



1 1 

1 2 

1 





10 
10 



-2 
-3 

1 



(7) 



(8) 



We use the inverses U 1 and V 1 to construct a unimodular transformation extending U 1 with 
the identity matrix, as follows: 



-3 2 
-2100 
10 
1 



1 1 

12 

1 





(9) 



and this product gives the transpose of M, the matrix in the unimodular transformation of (JH). 
This examples illustrates the case when U is not the identity matrix and where we may ignore S 
as its diagonal elements are all equal to one. 

We point out that the vectors in the null space of the exponent matrix A as in ([3]) are typically 
normalized so that the greatest common divisors of the components of the vectors equals one. 
We may change coordinates so that the first vector in the null space has only its first coordinate 
different from zero, the second vector in the null space can have nonzero entries only in the first 
two coordinates, etc. 
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Although we prefer to represent the solution set using only integer exponents for the param- 
eters, this is not always possible, consider for example 



B 



2 6 17 9 
4 14 13 3 



(10) 



The divisors for the two rows of B (and the denominators of the exponents of the parameters) are 
obtained via the Hermite normal form of B: UB = H, where U is a square unimodular matrix 
and H an upper triangular matrix. We assume that B is full rank and that the columns have 
been permuted so H has only nonzero elements on its diagonal. Let D be a diagonal matrix of 
the same dimensions as U which takes its elements from the corresponding diagonal elements of 
the matrix H. Then the coordinate transformation is defined by 



M 



D~ X B 
I 



(11) 



where / is the identity matrix, 
extended unimodular matrix 



To show that the determinant of M equals ±1, consider the 



U 



U 








(12) 



Because U is unimodular, det(f7) = ±1 and det(fJM) = ±det(M). We have det(UM) = ±1, 
because UM is an upper triangular matrix with ±1 on its diagonal as a result of the multiplication 
by D-\ 

Note that the rational exponents will appear only in the powers of the parameters as per- 
forming the coordinate transformation x = y M on the system x A — c = eliminates the first d 
variables of the ci-dimensional solution set. 



2.3 Solving Binomial Systems 

We denote a binomial system by x A — c = 0, where A G z kxn and c = (cq, c±, . . . , Ck-\) T with 
Cj for all i = 0, 1, . . . , k — 1. If the rank of A equals k, then k is the codimension of the 
solution set. Given the tuple (A, c), the solution set of x A — c = is described by a unimodular 
transformation M and a set of values for the last n — k variables. 

In the sketch of the solution method below we assume that A has rank k, otherwise x A — c = 
has no (n — /c)-dimensional solution set for general values of c. The steps are as follows: 

1. Compute the null space B of A, d = n — k. 

2. Compute the Smith normal form (U, S, V) of B. 

3. Depending on U and S do one of the following: 

• If U is the identity matrix, then M = V^ 1 and the first d variables have positive 
denominators in their powers when not all elements on the diagonal of S are equal to 
one. 
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• If U is not the identity matrix and if all elements on the diagonal of S are one, then 
extend U~ 1 with an identity matrix to obtain an n-by-n matrix E that has U~ l in its 
first d rows and columns. Then, M = EV . 

• In all other cases, define M as in (|lip . 

4. After the coordinate transformation x = y M , compute the leading coefficients solving a 
binomial system of k equations in k unknowns. Return M and the corresponding solutions 
of the binomial system. 

The solution procedure for binomial systems outlined above returns a representation with d 
parameters for the d-dimensional solution set which can geometrically interpreted as follows. For 
zero values of the parameters, we obtain the points of the solution set intersected with the first 
d coordinate hyperplanes. For nonzero values of the parameters, the powers of the parameters 
correspond to a choice of the basis for the null space of the exponent matrix of the binomial 
system. 

3 Sparse Polynomial Systems 

To look for d-dimensional components of sparse polynomial systems, we investigate solutions of 
initial forms defined by cones of normal vectors. In order for the initial form systems to have 
solutions with all coordinates different from zero, they need to be at least binomial systems. 

Although not all (and perhaps only few) initial form systems are binomial, the unimodular 
transformations explained in §2.2 are applied on a matrix of pretropisms. 

3.1 Initial Forms 

A polynomial / in n variables x = (xq, x\, . . . , x n -\) is denoted as 

/(x) = ^c a x a , c a €C\{0}, (13) 
aeA 

x a = Xq°x°^ ■ ■ ■ x^Si , where A is the set of all exponents of monomials with nonzero coefficient. 
The set A is the support of / and the convex hull of A is the Newton polytope P of /. Any nonzero 
vector v defines a face of P, spanned by 

m v (A) = { b G A | (b,v) = min(a, v) }, (14) 

agA 

where (•, •) denotes the usual inner product of two vectors. We use the notation in v (A) because 
a face of a support set defines an initial form of the polynomial /: 

in v (/)(x)= c a x a , (15) 

a6in v (A) 

where A is the support of /. For a system /(x) = and a nonzero vector v, the initial form 
system in v (/)(x) = is defined by the initial forms of the polynomials / with respect to v. 

Because the initial coefficients of Puiseux series expansions are solutions to initial form sys- 
tems, the initial forms we consider must have at least two monomials, otherwise the solutions will 
have coordinates equal to zero and are unfit as leading coefficients in a Puiseux series development. 
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3.2 An Illustrative Example 



In this section we indicate how the presence of a higher dimensional solution set manifests itself 
from the relative position of the Newton polytopes of the polynomials in the system. To illustrate 
a numerical irreducible decomposition of the solution set of a polynomial system, the following 
system was used in 



f(?,y,z) 



(y - x 2 )(x 2 + y 2 + z 2 - l)(x - 0.5) 
- x 3 )(x 2 + y 2 + z 2 - l)(y - 0.5) 








(16) 



, (y — x 2 )(z — x 3 )(x 2 - 
The solution set Z = / _1 (0) is decomposed as 

Z 



y 



z 2 - l)(z - 0.5) = 



Z 2 U Zi u Z Q 
{Z 21 } U {Z n U Z 12 U Z 13 U Z 14 } U {Z i} 



(17) 
(18) 



where 



1. Z 2 i is the sphere x 2 + y 2 + z 2 — 1 = 0, 

2. Z u is the line (x = 0.5, z = 0.5 3 ), 

3. Z\2 is the line (x = v0.5, y = 0.5), 

4. Z13 is the line (x = — yCU), y = 0.5), 



5. Z14 is the twisted cubic (y 



0,z 



x 



0). 



6. Zqi is the point (x = 0.5, ?/ = 0.5, z = 0.5). 



A first cascade of homotopies in [23] needed 197 solution paths to compute generic points on all 
components. The equation-by-equation solver of [27] reduced the number of paths down to 13. 
The Newton polytopes of the polynomials in the system are displayed in Figures [1] and (2) 

Consider a point on the 2-dimensional solution component of / _1 (0) and let the first coordi- 
nate of that point go to zero. As x\ = t — > 0: 



' m (i,o,o)(f)(x,y,z) 
( y(y 2 + z 2 
,2 1 ^2 



Alternatively, as x 2 = s 



1)(_0.5) = 
z(y 2 + z 2 - l)(y - 0.5) =0 
I yz{y 2 + z 2 - l)(z - 0.5) =0. 

-7- 0, we end up at a solution of the initial form system: 

'm( 0> i )0 )(f)(x,y,z) 

-x 2 (x 2 + z 2 -l)(x-0.5) =0 

-< (z-x 3 )(x 2 + z 2 -l)(-0.5) =0 

-x 2 (z - x 3 )(x 2 + z 2 - l)(z - 0.5) = 0. 



(19) 



(20) 
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Figure 1: From top to bottom, we see the Newton polytopes of f± and /2, the polynomials in (|16|) . 
The edges of the faces of the polytopes with normals (1,0,0) and (0,1,0) are marked in bold, 
respectively in red (thick solid lines) and black (thick dashed lines). 



Looking at the Newton polytopes along v = (1,0,0) and v = (0,1,0), we consider faces of 
the Newton polytopes, see Figures CD and [2j 

Combining the two degenerations, we arrive at the initial form system: 

in(o,i,o)(in(i,o,o)(/))(^2/^) 

y(z 2 -l)(-0.5) = 

_ l)(-0.5) = (21) 

yz(z 2 - l)(z - 0.5) = 

The factor z 2 — 1 is shared with inn i o,o)(i n (o,i,o)(/))( a; ) z )- 

Based on these degenerations, we arrive at the following representation for a solution surface. 
The sphere is two dimensional, x and y are free: 




1 + cotl + atl 



(22) 



For to = and t\ = 0, z = 1 is a solution of z 2 — 1 = 0. Substituting (x = to,y = t\,z = 
1 + cotg + c i^i) i n to the original system gives linear conditions on the coefficients of the second 
term: Co = —0.5 and c\ = —0.5. 

3.3 Asymptotics of Algebraic Surfaces and Puiseux Series 

Denoting by d the dimension of the algebraic surface defined by /(x) = 0, for x G C n , we assume 
the defining equations are in Noether position so we may specialize the first d coordinates to 
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Figure 2: The Newton polytopes of the third polynomial in (|16p . The edges of the faces of the 
polytopes with normals (1,0,0) and (0,1,0) are marked in bold, respectively in red (thick solid 
lines) and black (thick dashed lines). 



random complex numbers in /(x) = and obtain a system with isolated solutions. Moreover, we 
assume that when specializing the first d variables to zero, the algebraic set remains of dimension d. 
Geometrically this means that we assume that the algebraic set meets the first d coordinate planes 
(perpendicular to the first d coordinate axes) properly. 

We consider what happens when starting at a random point on the surface we move the first 
d coordinates to zero. For simplicity of notation we take d = 2 and consider a multiparameter 
family of polynomial systems: 



with co,ci,ci 5 i E C \ {0}, t>o,i; u i,i ^ Qj letting to and t\ go from 1 to 0, starting at a generic 
point on the surface with its first two coordinates equal to cq and c\. 

The multiparameter family in (|23p specifies the last equation as a series to leave enough 
freedom for the actual shape of the surface. While we may always move xq as going linearly to 
zero, with xq = coto, the second coordinate of a point along a path on the surface may no longer 
move linearly. Taking x\ as c\t\ would be too restrictive. 

As we move xq to zero as to goes to zero, then x\ can go to zero as well if «o,i > and 
Vi 7 i > 0, or go to infinity if vq,i < or vi t i < 0, or go to c\Ci t \ if both vq 7 i = and v\ t \ = 0. 
The multiparameter family in (j23H contains what we define as a multiparameter version of a 
Puiseux series for algebraic curves. Similar to x±, the other components of the moving point can 
be developed as a generalized Puiseux series 






Goto 




(23) 



x k = c k tl ' k t v 1 1 ' k {c 1>k + 0{toM)) 



(24) 
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Ck, ci,k G C \ {0}, VQ t k, vi : k G Q. If in the limit — when to and t\ are both zero — the solution is 
finite and of multiplicity one, and if the powers in the series are positive integer numbers, then 
the generalized Puiseux series coincides with a multivariate Taylor series. 

As to and t\ go to zero, the system F(to,ti) = — obtained after replacing xq and x\ using 
the last two equations of (123f) and after substituting (|24jl for the remaining n — 2 into /(x) = — 
must have at least two monomials with lowest power in to and lowest power in ti in every equation 
because c^, G C \ {0} for all k = 0, 1, . . . ,n — 1. We call the part of /(x) = corresponding 
to F(to,ti) with lowest powers of to and t\ the initial form system of /(x) = with respect to 
the normal vectors vo = (l,«o,l> Wo,i, . . . ,«o,n-i) and vi = (0, «x,l> w i,2, • • • , Wi in _i). Because the 
normal vectors are the leading powers of the generalized Puiseux series, vo and vi can be called 
tropisms in analogy to the case of algebraic curves. 

The derivation of Puiseux series for an algebraic set in any dimension d if formulated as 
follows. 

Proposition 3.1. If /(x) = is in Noether position and defines a d-dimensional solution set 
in C n , intersecting the first d coordinate planes in regular isolated points, then there are d linearly 
independent tropisms Vo,Vi, .. .v^-i G Q n so that the initial form system in vo (in vl (- •• in Vd _ 1 (/) •• • 
y M ) = has a solution c G (C \ {0}) n . This solution and the tropisms are the leading coeffi- 
cients and powers of a generalized Puiseux series expansion for the algebraic set: 

x = t 

X\ — Tq i x 

j^O.d-l .Vl,d-1 ,Vd-l,d-l 

Xd-1 — l z l ' ' ' l d-l 

(25) 

Xd = c t t x ■■■ t d _ 1 H 

a^d+i — cit t x • • • t d _ 1 H 

.«0,n-l .Vl,n-1 ,fd-l,n-l , 

Xn — C n —d—\i'o ' d— 1 ■+"""" 



Proof. Because the set defined by /(x) = is in Noether position, we can let the first d variables 
go to zero, using for example a multiparameter homotopy as in (|23p and still obtain regular 
isolated solutions, denoted as (0, 0, . . ., 0, cq, ci, . . . , c n -d-\) G C n . 

The tropisms Vq, vi, . . ., Vd-i define the initial form system, i.e.: those monomials in the 
system /(x) = that become dominant as the parameters to, ti, t^_i move to zero. In 
particular: for any vector v in the cone spanned by the tropisms, we have that every monomial 
x a in the initial form system makes minimal inner product (a, v), minimal with respect to any 
other monomial x b not in the initial form system, i.e.: (a, v) j (b,v). 

Because the leading terms of the Puiseux series vanish at the initial form system, the inner 
product with the monomials and the leading powers must be minimal compared to all other 
monomials in the system. Hence the shape of the Puiseux series. □ 
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3.4 Polyhedral Methods 



In our algorithm to develop Puiseux series developments for algebraic sets, Proposition 13.11 is 
applied as follows. If we are looking for an algebraic set of dimension d and 

• if there are no cones of vectors perpendicular to edges of the Newton polytopes of /(x) = 
of dimension d, then the system /(x) = has no solution set of dimension d that intersects 
the first d coordinate planes properly; otherwise 

• if a d-dimensional cone of vectors perpendicular to edges of the Newton polytopes exists, 
then that cone defines a part of the tropical prevariety. 

We call a vector perpendicular to at least one edge of every Newton polytope of /(x) = a 
candidate tropism or pretropism. 

Algorithms to compute a tropical prevariety are described in [7]. As we outlined in JT], we 
applied cddlib [13] to the Cayley embedding of the Newton polytopes of the system to compute 
pretropisms. With the Cayley embedding we managed to compute all pretropisms of the cyclic 
12-roots problem, reported in pQ. 

For highly structured problems such as the cyclic n-roots problem, a tropism found at lower 
dimension often occurs also in extended form for higher dimensions. For example, for n = 4, a 
tropism is (+1, —1, +1, —1) which extends directly to (+1, —1, +1, —1, +1, —1, +1, —1) for n = 8 
and (+1, — 1, +1, — 1, +1, — 1, +1, — 1, +1, — 1, +1, — 1) for n = 12, and any n that is a multiple 
of 4. 

In addition to the extraneous results reported from the Cayley embedding, it suffices to restrict 
to pretropisms with positive first coordinate because geometrically we intersect the solution set 
with the coordinate hyperplane perpendicular to the xo-axes at the end of moving xo to zero. 
Allowing a negative first exponent in the first pretropism corresponds to intersecting the solution 
set at infinity, when in the limit we let xq go to infinity. 

In any case, after the computation of pretropisms, exploiting permutation symmetry is rela- 
tively straightforward as we can group the pretropisms in orbits and process only one generator 
per orbit. 

3.5 Puiseux Series for Algebraic Sets 

The approach to develop Puiseux series proceeds as follows. For every d-dimensional cone C of 
pretropisms: 

1. We select d linearly independent generators to form the d-by-n matrix A and the corre- 
sponding unimodular transformation x = y M . 

2. Because the matrix A contains pretropisms, the initial form system in vo (in Vl (• • • in Vd _ 1 (/) • • • ))(x) = 
determined by the rows vo, vi, . . ., v^-i of A has at least two monomials in every equa- 
tion. If the initial form system has no solution with all coordinates different from zero, then 

we move to the next cone C and return to step 1, else we continue with the next step. 
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3. Solutions of the initial form system found in the previous step may be leading coefficients 
in a potential Puiseux series with corresponding leading powers equal to the pretropisms. If 
the leading term satisfies the entire polynomial system, then we report an explicit solution 
of the system and we continue processing the next cone C. Otherwise, we take the current 
leading term to the next step. 

4. If there is a second term in the Puiseux series, then we have computed an initial development 
for an algebraic set and report this development in the output. 

To compute in the last step a second term in a multivariate Puiseux series seems very com- 
plicated, but we point out that it is not necessary to compute the second term in all d variables. 
To ensure that a solution of an initial form system is not isolated, it suffices that we can compute 
a series development for a curve starting at that solution. In practice this means that we may 
restrict all but one free variable in the series development and apply the methods we outlined 
in pQ for the computation of the second term of the Puiseux series for a space curve. 

With Puiseux series, the solutions of the initial form system can be extended to form a 
witness set. A witness set |26j is a numerical data structure for positive dimensional solution sets 
of polynomial systems. Depending on the heights of the powers in the series, we may need more 
than the second term to ensure convergence with Newton's method. 

4 Applications 

Our polyhedral approach enables to compute exact representations for positive dimensional so- 
lution sets of the cyclic n-roots problem (pQ). 

4.1 On cyclic 9-roots 

Taking n = 9 in ([I]), for cyclic 9-roots, we show that our solution can be transformed into the 
same format as in the proof we found in [12\ Lemma 1.1] of the statement in [2] that square 
divisors of n lead to infinitely many cyclic n-roots. 

Among the tropisms computed by cddlib [13] on the Cayley embedding of the Newton 
polytopes of the system, there is a two dimensional cone of normal vectors spanned by u = 
(1, 1, —2, 1, 1, —2, 1, 1, —2) and v = (0, 1, —1, 0, 1, —1, 0, 1, —1). The vectors u and v are tropisms. 
The initial form system in u (in v (/))(x) = is 
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X 2 + X 5 + X 8 = 

X Xg + X 2 X 3 + X 5 X 6 = 
XqX\X2 + XqXiXs + X X 7 Xg + X1X2X3 
+ X 2 X 3 X 4 + X3X4X5 + X 4 X 5 X 6 + X 5 X 6 X 7 

+ X 6 X 7 X 8 = 

X XiX 2 X 8 + X2X3X4X5 + X 5 X 6 X 7 X 8 = 

X XiX 2 X 3 Xs + XoX 5 X 6 X 7 X8 + £2X3X4X5X6 = 

X0X1X2X3X4X5 + X0X1X2X3X4X8 

< + X0X1X2X3X7X8 + XoXiX2X 6 X 7 X8 (26) 

+ XoXiX 5 X 6 X 7 X8 + XoX4X5X 6 X 7 X8 
+ XiX2X3X4X 5 X 6 + X2X3X4X 5 X 6 X 7 

+ X3X4X5X6X 7 X8 = 
X0X1X2X3X4X5X8 + XoXiX2X 5 X 6 X 7 X8 

+ X2X3X4X 5 X 6 X 7 X8 = 
XoXiX2X 3 X4X 5 X 6 X8 + XoXiX2X3X 5 X 6 X 7 X8 

+ XoX2X3X4X 5 X 6 X 7 X8 = 
XoXiX2X3X4X5X6X 7 X8 — 1 =0. 

Although not binomial, in u (in v (/))(x) = is is significantly sparser and thus easier to solve than 
the original system. To solve in u (in v (/))(x) = 0, we eliminate xo and xi with a unimodular 
coordinate transformation M that has u and v on its first two rows. The last seven rows of M 
are zero except for the ones on the diagonal: 



" 1 


1 


-2 


1 


1 


-2 


1 


1 


-2 





1 


-1 





1 


-1 





1 


-1 








1 





























1 





























1 





























1 





























1 





























1 





























1 



(27) 



The matrix M defines the unimodular coordinate transformation x = y : 

xo = Vo x 3 = y y 3 x 6 = y y 6 

xi = yoyi x 4 = 2/oyiy4 x 7 = y yiy 7 (28) 

X2 = y 2 Vi 1 y2 x 5 = y^y^ys x 8 = y^y^ys- 

The transformation x = y M reduces the initial form system in u (in v (/))(x = y ) = to a system 
of 9 equations in 7 unknowns. 

After adding two slack variables to square the system (see [26J for an illustration of introducing 
slack variables), the mixed volume equals 326. In contrast, the mixed volume of the original 
polynomial system equals 20,376. 

We find that the entire cyclic 9-roots system vanishes at this first term of the series expansion. 
Recognizing the numerical roots as primitive roots of unity leads to an exact representation of 
the two dimensional set of cyclic 9-roots. 
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Denoting by u = e* 27r / 3 the primitive third root of unity, u 3 — 1 = 0, our representation of the 
solution set is 

x = t x 3 = t u Xq = t u 2 

x\ = *o*i X4 = totiu X7 = totiu 2 (29) 

Introducing new variables yo = *0) Hi = *o*i> and 1/2 = *q 2 t^ 1 u 2 , our representation becomes 

xo = yo x 3 = y u x 6 = y u 2 

x\ = yi x 4 = yiu x 7 = y x u 2 (30) 
x-i = 2/2 x 5 = y 2 u x 8 = y2U 2 

which modulo y^yfy^u 9 — 1 = satisfies by plain substitution the cyclic 9-roots system, as in the 
proof of [12 s Lemma 1.1]. 

Note that the representation in (|29p allows a quick computation of the degree of the surface. 
This degree equals the number of points in the intersection of the surface with two random 
hyperplanes. Using (|29p for points on the surface, the two random hyperplanes become a system 
in the monomials to, io*ii an d *o ^1 '■ 



2 1 

aito + ai,2*o*i + ci-2,-i*o h =0 (oi\ 
Pito + /hfltoh + /3_2,-l* V = 

for some complex numbers 014 j and The above system is equivalent to the system 

*o *i - co = , > 

*i - Cl = [62) 

for some Co, c\ G C. We see that for any nonzero Co and ci, the system has three solutions. So the 
algebraic surface represented in ([29]) is a cubic surface. Using other roots of unity and permuting 
variables leads to an entire orbit of cubic surfaces. 

Using the representation (|30p . we arrange the position of the coefficients with it as a third 
root of unity: 



1 


u 


u 2 


u 


u 2 


1 


u 2 


1 


u 


u 2 


u 


1 


u 


1 


u 2 


1 


u 2 


u 



(33) 



shifting the variables in forward and backward order. So the one cubic surface leads to an orbit 
of 6 cubic surfaces, corresponding with our numerical results of [25]. 

4.2 On cyclic m 2 -roots 

While the Cayley embedding becomes too wasteful to extend the computation of all candidate 
tropisms beyond n = 12, by the structure of the tropisms for n = 9 we can predict the tropisms 
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for cyclic 16-roots: 



u = (1, 1, 1, -3, 1, 1, 1, -3, 1, 1, 1, -3, 1, 1, 1, -3), 
v = (0, 1, 1, -2, 0, 1, 1, -2, 0, 1, 1, -2, 0, 1, 1, -2), 
w = (0,0, 1,-1, 0,0, 1,-1, 0,0, 1,-1, 0,0, 1,-1), 



(34) 



and the corresponding initial form solutions are primitive fourth roots of unity. Similar to (I29p 
and (|30p we can show that the exact representation obtained with tropical methods corresponds 
to what is in the proof of \12\ Lemma 1.1]. 

A general pattern for surfaces of cyclic m 2 -roots is below. 

Proposition 4.1. For n = in 2 , there is an {m — 1)- dimensional set of cyclic n-roots, represented 
exactly as 

x km+0 = V>kto 
Xkm+1 = Uktotl 

Xkm+2 = Uktotlt2 , „ 

(35) 



for k = 0, 1, 2, . . . , m — 1 and Uk = e 



■£fcm+m— 2 
%krn+m— 1 — 'Ufc^o 
i2kir/m 



Uktot\t2 ■ ■ ■ t m -2 



m+1 j.— m+2 
l l 



l m-Z l m-2 



The substitution to = so, *o*i = s±, £0*1^2 = S2, • • •, t 
simplifies ([35]) . 



m+l i— m+2 



■ ■■r 2 t- 1 



-1-1 



S m-2 



Proposition 4.2. Tfte (m — 1)- dimensional solution set in (|35p /ias degree equal to m. 

Proof. To determine the degree of an (m — l)-dimensional algebraic set, we intersect the set with 
m — 1 hyperplanes with random coefficients. In any linear equation we replace the x-variables 
using the equations in ([35]) . dividing each equation by to to obtain a nonzero constant coefficient. 
Because every Xj corresponds to one monomial in to, t\, . . ., t m —i, bringing the coefficient matrix 
into a reduced row echelon form leads to a binomial system of m— 1 equations inm-1 unknowns: 



m + — m+2^— m+3 



t 



■m-2 C 
h ~ Cl 

ht2 - c 2 











(36) 



M2*3 - - - t m -2 — C m -2 

Collecting the coefficients {cq,c\,C2, ■ ■ ■ , c m -2) in c and the exponents in a matrix A, we denote 
the binomial system as t A = c with 



A 



-m —m + 2 — m + 3 












1 

1 
1 



-2 -1 





1 

1 1 



(37) 
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The binomial system has | det(^4)| = m solutions and therefore the degree equals m. □ 

Applying the permutation symmetry, shifting the variables forward and backward as in (j33|) . 
we find 2m components of degree m. 
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